Approved For Release ST AT 



Approved For Release 


2009/08/26 : 

CIA-RDP88-00904R0001 001 1 0 




I 


Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 0011 001 2-6 


i 


w 


Third United Nations 
International Conference 
on the Peaceful Uses 
of Atomic Energy 


A/CONF.28/P/364 

UUSH 

May 1964 

0:-iiiinal: IJUSSIAI.’ 


Confidential until official release during Conference 


SOME m-l ASPECTS OF THE APPLICATION OF THE ADJOINT 
FUNCTIONS AND OP THE PERTUPDATION TirEORY IN 
REACTOR AND SHIELDING DESIGN. 

A.A.Abagian, G. I. DruJ Inina, A. A. Dubinin, S. W. Zarltsl^y , V.V. Orlov, 
V.Ya.Pupko, A. P. Suvorov, L.N.Usachev, R. P. Fedorenko. 

The application of perturbation theory formulas for a critioal nuc- 
lear reactor, published e.g. in Refs /1/, /?/, /3/, /A/, /5/, proved to 
be rather suitable in reactor deoi£pi and experimental research. 

Refs /6/, /?/ extended the area of perturbation theory application 
for arbitrary linear lunotionals in stationary and time depending prob- 
lems of the radiation transfer with given souroes. In Ref /8/ the per- 
turbation theory formula is obtained for a ratio of different processes 
numbers in a stationary reactor. 

The application of formulas from Ref/ 8/ giving more possibilities 
in an investigation of reactor physical characteristics is discussed in 
chapter I (L.N.Usachev, S.M. Zaritsky). 

The ezampl'' of this theoiy applying for breeding ratio calculations 
for variants of fast breeder is presented. Th.e possibility of using this 
theory for improvement of the multigroup constants system on the basis 
of results measured in a critioal assembly is considered. 

The use of perturbation theory formulas for the effective multipli- 
oatlon constant and for latio of different processes numbers in a reactor 
for the statement of variational problems on optimum distribution of dif- 
ferent materials in a reactor with an arbitrary neutron spectrum is dis- 
cussed in chapter II. (V.Ya.Pupko, G. I.DruJ inlna). Some identical prob- 
lems for in applying to a thennal reactor were considered in Ref/^*/. 

Chapter III (V.V. Orlov, A.A.Abagian, R. P. Fedorenko, A. A. Dubinin, 

A. P. Suvorov) is devoted to using of perturbation theoiy formulas for 
the shielding characteristics optimization. 

As the application of olassicsl oaloulus of variations methods for 
shielding problems proved to be difficult in some cases, the optimiza- 
tion process is made on the base of linear programming ideas. 

[25 YEAR RE-REVIEWl 


Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 001 10012-6 





Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 0011 001 2-6 


Gh apter 1» THE PERTUHBATION THEOW EuR A ILM’TO OF 

PROCESSED RUMBEUS IN A REACTOR. 

Duimnary and some improvementB in the derivation 
made in Ref /8/. 

The well known in the reactor theory perturbation theory it 'Dlion 
for determination of disturbed critical reactor parameters (see for 
instance /3/) 

1 - J(- = c:/£cJfy}/ ^ ^ ^ 

is obtained from the equation for neutron flux ?*(?, E, "n) in a distur- 
bed reactor with primed parameters: 

and from the equation for neuti'on importance in an undisturbed reactor 


-Dfj. yf'F.l: '''! / IF'ffe'n'j 

To obtrviu the relation (1) it is neoeBeary to cross multiply equa- 
tions / 2 / and /3/ by and P*, cubtraot one from another and integrate 
over all energies, velocity directions and reactor volume up to the 
outer boundary with vacuum, 

Identloal operations with co-adjoint equations, formulated in Rejf 
A/> were used in Ref /?/ to obtain the perturbation theory relation 
for an arbitrary linear functional of neutron flux. There was conside- 
red the more general time depending problem on the distribution of neut- 
rons from given sources. 

The concept of the neutron importance in any previous cycle to 
a process in a chosen (zero) oycle was introduced in Ref /8/ to ob- 

tain the perturbation theory relation for a ratio of processes numbers 
ocourred in a stationary reactor which is described by homogeneous equ- 
ations. The set of equations for these neutron import anoes serves as 
adjoint equation in derivation of the relation required. If we are in- 
terested, for example, in the neutron capture in i^^ isotope this set 

O' ->^'-'-'>17?,. 






Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 001 10012-6 




Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 0011 001 2-6 


Th6 BuooesBlve oolving of equations (4a), (4b), (4o) is equivalent 
to a usual procedure of solving of equation (3) by the suooessire Itera- 
tive method. Henoe, it is oloar, that at A Using the ad- 

JointneSB of left parts of equations (3), (4a), (^b)^ (4o) and left 
part of the equation for the undisturbed flux P(ri E, “n) (equation /2/ 
being without primes) and carry ing out the procedure of the cross multi- 
plying and subtraction, as It was indloated above, v/e obtain the set of 
equations; 

Qi .... (4d) 

K~-7~ » where ^ 


Henoe 


for 


Summarizing equations for P'tf Pt*.# to P^, we obtain the equation 


Carrying out the same procedure of the oross multiplying, integra- 
tion with equations C2) and C5), then forming the ohunge of the number 
of prooesses under consideration which is equal to 
cTf?,. = a/ -Qi etSJSlc^V 

and using the relation C^d): 

c^FcISl^U = P/F0/9 

we obtain ^ ✓ /- j* / 

- * /fF.^Cf’-F)JO ( 6 ) 

^Considering only pervurbatlons which remain reactor oritloal, l.e. 
perturbations satisfying the condition Cl), we note that in oeoond, 
third and fourth terms of the right part of (6) it ie porrlble to re- 


place which tends to infinity at n by 

which has the definite limit 










(7) 


The last term in the right part is equal to 

'S yF'- Fjc^G 

Only ratios of prooesses numbers (for instance the ratio of oapture- 
to-fission numbers, i.e. to the reactor power) as well as changes of 
these ratios when making perturbations in reactor have the physical mea- 
ning in a critical reactor. 

The variation of ratio of capture number in Isotope- 1 to this in 
Isotope-j which equals 


■Oi I Gc _ Oc _ Q,- /'(Ta.- _ 

" d- ~ oTt a a,- J 


mi- 




"j I ^ ^ y 

iB possible to write using (s) at n^oo and(7) in the following fora 

364 . 3 _ 
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^(u)'- vye - /(^ -'Mjf'Y' ( 9 ) 

Here 

Is the neutron flux in a disturbed reactor normalized 
in such a way that capture number in isotope-J is equal to this number 
in an xmdisturbed reactor, and 

^ Qi ay 

is a function having the meaning of neutron Importanoe with respect to 
ratio of processes numbers as it is seen from relation (^9)* 

Prom relation (11) and the linearity of set^of ^quations (zi) It is 
clear that eubstitutlng the linear oomhlnatlon^^;^ Instead of 

^/Fa in the right part of equation (lia) and Buooesalvely solving 
equations (h) we shall directly obtain a sum ^ . Only one 

set of equations is solved instead of two and there is no need for sub- 
tracting the main harmonic P^ from each as it is done in relation (7 ), 

since it is generated by every item in the right part of the equation 
with the same coefficient 1/1^ but with a different sign. This combina- 
tion which is advisable to use in oaloulatlons was pointed out In Ref 13 
The equation for the neutron importance with respect to linear-fractional 
functional was derived In this reference based upon the introduced 

there neutron importance oonoept with reepeot to an arbitrary non- 
linear functional and upon the neutron importanoe balance wbioh had 
been used in Ref 3 and ?• Por small disturbances it was shown 
that the functional derivative of the mentioned non-linear functio- 
nal with respect to neutron flux should be in the right hand of the 
equation for the neutron importanoe. In the case of the linear- 
fraotional functional, under oonsldfiratlon, this derivative is equal 

The algorithm of the perturbation theory described was put in nrao- 
tloe by us using 18-group diffusion approximation in the program for 
digital computer. All examples of the perturbation theory application 
^^^esorlbed below have been realized or are being realized now by means 
of this program. 

§2. The examples of the use of perturbation 
theory formulas. 


A a Reaotor variant oaloulations. Reactor design variations are 
usually described, from the neutron physios standpoint, by change of 
the material oonoentratlon in reaotor different points. Henoe, for reactor 
variant oaloulations by means of formulae (l) and ( 2 ) reaotor parameters 

- 4 - 
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^ - 21 S'P a; 


-_iL 

variations are presented as: ^a4ys ^s/s 

where dfefine the Interaction between neutrons and '^-type 
nuoleus, / 

To show -that perturbation theory formulas are rather suitable in 
the applying Tor variant oalculations we shall write them for a given 
reactor. Let us oonsider a fast reactor which is the spherloal simpli- 
fied model of a power breeder. In the core with a radius 66 cm and with 
a volume 1204 litres fuel, i.e. uranium and plutonium oxide, occupies 
35* of all volume, steel - 25* and sodium - 40*. Nuclear concentrations 
of oxygen, plutonium, uranium, steel and sodium (multiplied by in 

1 om^ are correspondingly equal to ^*^=0.01288, *0,001 100, 

0.005334, b0.0214, *0,00888. 

In the breedi^ blanket with a thickness 60 cm concentrations are: 
-0.0272, ^-lo, -0.01355, ^"'-0.0153. -0.00/153. 

In ■^.e steel reflector with a thickness 20 om. -0, ^-'-0, 

P- =0.072. P- =0.00345. The total breeding ratio BR and the breeding 

<JS$ ^ 

Of the core BRC are correspondingly equal to 1,4575 and 0.5947. 

At unifonn variations of different materials concentrations in the 
core and in tbe breeding blanket it is suitable to present relations 
describing variations of the plutoniiim oritioal ooncentratlon,BH and 
BRC in the next form: ooon 

here N Is* the reactor zone number. Coefficients formulas 

(12^ are presented in the table: 


K 

p^239 y238 

Na 

Steel 

0 


c 

Fission 

fraRments 

Mo 

Hi 

- 


0.0411 



-0.0^132 

-0.o\41 

-0. 

0^194 

-0.0^191 

0.0689 

0,0204 

-0.0^168 

- 

0.119 

0.0^56 

-0.0^104 

-0.0^137 

-0. 

0^795 

-0.0^480 

-0.0^674- 

•O.O^IoI 

>-O.0''223 


884 

68.1 

-9.44 

-8.15 

-8. 

73 

-10.6 

-7.95 - 

•8,84 

-11,5 


303 

5.00 

-1.76 

-3.25 

-0. 

948 

-0.562 

-53.7 - 

• 17.6 

-5,44 


978 

102.4 

5.45 

2.27 

5. 

22 

6.95 

-0.767 

1.05 

1,60 

%C^ 

..99 

-0.363 

0.284 

0.214 

0. 

255 

0.281 

0.382 

0.163 

0,422 


These coefficients were calculated for small disturbanoes, i.e. 
in assumption ?»«P and P**«P. The mentioned progrannne permits to make 
the calculation with F» and P'* to determine llmitB of these formulas 
application and to clear up causes of possible disagreements. 

364 - 5 - 
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D. Effect of the constant 3 inaccuracy. At the calculations of ooef- 
floients for relations (l?)we obtain also ooeffloients defining the ef- 
fect of any material constant inaccuracy in any group on the oritioal 
concentration, total breeding ratio and core breeding ratio. It is pos- 
sible to perform the analysis similar to that performed in Ref l\k/ by 
means of these ooeff lolents. It will be recalled that mentioned coeffi- 
cients in Hef /14/were determined by the direct calculation. 

C. Analysis of critical assembly experiments for the constant system 
improvement . Perturbation theory equations (I) and (^ 9 ) may be con- 
sidered as linear equations, relating differences of observed and calcu- 
lated values with inaccuracy in the multlgroup constant system, which, 
as it is proposed, causes the calculation and the experiment dlsorepancy. 
The application of the perturbation theory relations with above meaning 
for analysis of oiitical assembly experiments and for suoh an experiment 
planning oan give the systematio improvement method of the, multlgroup 
constant system for a reaotor design. For each critical assembly on 
which experiments are performed it is possible to calculate the linear 
equation coefficients, relating the oritioal concentration variation and 
the changes of ratios of different prooeases numbers observed in an experi- 
ment to multlgroup con:-tajit variations. Having a large number of the 
mentioned relations for several critical assemblies and for many ratios 
of processes numbers it is possible to state the matching problem of 
constants, describing all experiments in the best way. In the general 
statement suoh problems oan be solved ly the linear programming methods 
/11/ and in particular simple cases they can be solved by the trial and 
error method. 

Capture-to-flssion number ratios in different Isotopes measured in 
different points and integrally over reaotor zones oan b <5 used as ratios 
of processes numbers. Suoh easily measured values as the reactivity at 
the placing of different samples oan be used too. Indeed, the reaotivl- 
ty oaused by a sample placed in a oritioal reaotor oan be >rritten so 
(so./ 3/) 

^ , .her. ( 13 ) 

is an additional neutron Importance arisen per u^it of time at the 
sample placing, 

is the importance of fission neutrons emlssioned throughout a reactor 
per unit of time. Hence, it is clear that to calculate reactivity va- 
riations ^ ) it is possible to use the same algorithm as one 
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for other ratios of prooeeses numbers. It is necessary to note that in 
this case values which at multiplying by neutron flux F and integrating 
over energies, velocity directions and reactor volume will give and 
ay determined from equations (13) and (ia), should be in right hand of 
equation (4a) instead of 1/la^ and l/ls^. Therefore, iv is neoessary to 
replace in equation (4a) and in following equations: 

Z/' f 

Formulas ( 9 ), (l5) and (l6}take into aooount not only possible In- 

accuracy in sample oonstants but inaoouraoies of the constant s describ- 
ing the x*eactor and also the reactor speotrum variations oaused by 
influence of the inserted sample. 

The variation of the prompt neutron mean life time being a ratio of 
the total reactor neutron importance to the Importanoe of fission neut- 
rons, emitted per second, oan be obtained Identical ly/j/. In this case 
in fonnula8(4)- (9j it is neoessary to replace# 
dEdnc/y ; Ve/ 

a; fpfj ^ d0 ; Vei-* J fJ ^ de dSL 

Such analysis^ was already begun tt) carry out at the treatment of 

experimental data obtained in aseecblleo |l5*/ 

The development of the perturbation theory desoribed ooneidorably 
extends its appliostion area which of course Is not limited hy presented 
examples (See e.g. /s/ ). In particular, it heoomes possible to set the 
problem of optimization of breeding ratio or some other ratios of proce- 
sses numbers which is discussed in the next chapter. 

C h a p t c r II. THB USHTO 07 IHS PSHTTOMTION THBORT 

PUHCTIOKAIS IB THE 7AHIATI0NAL PHOBUSMS 
OH HATEEIAL OPTHTOM DISTHIBOTION IN 
HUCLEAH HBACTOHS. 

Using the perturbation theory it is possible to state In general 
oase the problem on a different reactor materials optlJsua distribution 
with relation to the achievement of most advantageous oharaot eristics 
oaused by the effect of these nmterlals. The problem on reactor orltloal 
mass minimization and the problem on the maximum breeding ratio in a 
breeder are examples of suoh problems. The approach to the variational 
problem ooneists in the finding of extremum of corresponding funotlonal 
which has variation presented with perturbation theoiy formulae. Seve- 
ral other functionals or nonlntegrel conditions oan be simultaneously 
fixed so as from oaloulus of variations point of view we shall deal 
with the isoporimetrloal or conditional extremum problem /16/. In the 
Jilt - 7 - 
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first from mentioned problems the functional is the reactor effective 
multiplication constant having the variation being written with the 

perturbation theory formulas/ 3/. In the second problem perturbation 
theory formulas are used /3/ for the linear fractional functional which 
desorlbes a ratio of two different reaotor processes. Perturbation theo- 
ly formulas contain different materials effioienoy functions desorl^ 
blng the effect of these materials content variation on a value of vary- 
ing functional, This permits to fonnulate the variation problem In the 
optimum distribution investigation of different materials in a reaotor. 

The minimization reaotor charge problem and more briefly the maxi- 
mum breeding ratio problem are oonsidered below. Similar to considered 
problems it is possible to formulate other variation problems cn the op- 
timization of reaotor oharaoteristlos (the maximum oompensating oapaoity 
of the absorption material, eto,), 

§1, The statement of the pi^blem on minimum oritioal 
reactor mass. 

We have the nuclear reaotor with limited sizes and with the given 
content of structural materials. At some uniform oonoentratlon of U-235 
in the core the reactor has definite value of and oritioal mass 

G^. How should the distribution of U-235 over all core be modified to 
obtain the critical mass reducing without K ^ variations, if it Is po- 
ssible? 

Similar problem has been considered by Goetzel /l7/to the thermal 
reactor case with the same moderator in the oore as in the reflector. 

The method of suoh reaotor oaloulation is oonsidered also in monographs 
/13/ and /4/. 

In this report the oondltion of the oritioal mass minimization is 
oonsidered in general case for a reactor with an arbitrary neutron speo- 
truffl. 

To definltlbn we assume that varying the oonoentrption of U-235 ps’(^) 
in the point T we slmultaueouely change the content of only one compo- 
nent having the oonoentratlony? In general we have the possi- 

bllity In prlnolple to lower the oritioal mass of the given reaotor from 
^ 5 " 0 onst with the oorresponding redistribution of U-235 In the 
oore. The problem is oonsidered in two stages; at first the variational 
problem for maximum at the given reaotor mass Is stated and then 
the transition to problem on minimum G^ at fixed is disoussed. 

§2, The achievement of maximum ^ reactor with 

the given charge, 

YaryingJ) ^(r) we shall find the 

- 8 - 
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( 1 ) 


( 2 ) 


perimetrical conditions 

According to Euler's method /16/ we construct the functional 

^ (Ps) - ^•ff 

and find its absolute extremum. The expression for the variation is 
written so: 

The first term we shall express witfi perturbation theory formulas 

By symbol fC?| p ^ it ?s defined the net neutron Importtjioe generated 
per second in the point r with one fissionable material nucleus as a 
result of all interaction processes of neutrons with U-235. The function 
f(r5_^ c) la the following tart is named as effiolenoy function of U-235. 

DependLoe of this function upon P JX) is functional, its structure is: 

,n) P . (»> , 




where p - P v- " material oonblift rat ions corresponding to Its maxi- 
y J JC |0 

mum density . In ceneral: 

-Elist')] F dE'cln'dEc/n-l[F(^.£. si}FKBp)6a\% (e) dEdsi} 


*W '*’{£• 


/ '"’rvyo ) I ^ y/// F(7. e; H') h / F *a.E. A) dE'dndE cIR 




//* f- dK j 


wnere K ’is the normallzaf fbn importance or xission neutrons, 
fn 


( 6 ) 

IJuler^ s 


(7) 


equation for the variational problem under oonslderation 
is received by using equations (3) and ( 4 ) from ertremuin condition 

SJiPf] = 0 

Equation (7) is the necessary condition of an optimal distribution 
of the given amount of n-235 in the ooro, ensuring the maxlnrum in 

a reaotor with an arbitrary neutron spectrum . The physical meaning of 
this condition is clear: if the U-235 effiolenoy function is constant 
throughout the ooi*e, every gram of the flsBlonabla material in any point 
is used equivalently. Prom this it follows that hero we consider suoh reac- 
tors, for which an extremum of funotlonal exists at the oonoen- 

tration of U-235 O^/i 

Equation (7) permits to base the method of experimental detormlna- 


^ The same requirement for the minimization of the oritloal mase of a 
thermal reactor is noted in Ref /4| by A. Weinberg and E.Wigner. 

- 9 - 
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tlon of the U -235 ooncen 1 :ration distribution in i’eaotor for which maxi- 
nmm is realized at "the given charge. This equation can serve too as 

a proper solution criterion and should be used In numerical calculations 
for such reactors. These are accomplished with an iterative method, ac- 
hieving on prooedu’^3 reduction to constant value of efficiency function 
oX U-235 In a core volume. 

The possible oaloulating iterative procedure consists in gradual 
redistribution of U-235 in core in accordance with the shape of the oal- 
oulated efficiency function of U-235. 

§3* The transition from maximum problem to minimum 

charge problem. 

Let fimctlons ju ^(r) , p providing maximum 

have been found for different charge values of the reactor 

tmder consideration. Prom Pig, I it follows the reactor charge correspon- 
ding to curve 1 with ( ^ ^ 5)73 3 ^ ^ value of requi- 

red will be the least. In connection with this it is necessary specially 
to Analyse the sign of derivative ^ 3 different 

reactors. 

Using the perturbation theory formulas it oan be shown that: 

' d 'Va.l. VflLj J y 

The critical reactor mass minimization at / < 0 requires an additi- 
onal consideration. Identically to a proof in Ref /iO/ for reactors with 
p ^aconst it oan be, shown that for the different types of reactors with 
variable concentration of U-235 the details of function diagram 0^=* 
fCVa, 3 ) are determined with a relation the neutron leakage factor from 
a reactor and an efficiency fissionable material factor ? — f in a reactor, 
Schematioal diagrams of dependence G^af(Va,3) for the different spectrum 
reactors are given in Pig. 2, In all oases but the part CD of the curve 
the U-235 effiolenoy factor ^ >0, In the part CD a depreciation of 
U— 235( < 0 ) ooours and the slope of the curve G^»f(Vr.,3) varies. Causes 
of Buoh event oan be, for instance, a higher radiative capture of neut- 
rons in the fissionable material, an availability of the strong resonance 
absorbers In reeotor, the using of a hydrogenous moderator, eto. The curve 
slope change of G^»f(Va. 3) In part CD leads to that without varying the 
core volume certain amount of U-235 oan be removed from suoh reactors 
replacing it by a moderator without reactivity loss. Thus the critical 
reaotor mass minimization at ^^0 is reduced pro forma to transition 
from part CD to part C*D of the diagram curve, 

36^ 

- 10 - 
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§/u On fission density non-uniformity In minimized 
oharge reaotorn. 

For minimizing a charge It serves to purpose to displace U-235 
from peripheral reactor regions near to central regions, however such 
fission material redistributing in reactor is able adversely to effect 
on the non-uniformity coefficient of heat generation in a core. At cri- 
tical mass minimizing in a small reaotor it is neoesaary to raaKe a*^'^ooal' 
rifling of U-235 concentration on periphery of a core (Ref 17) parallel 
with the fissionable material redistributing to reactor center and some 
improvement of fission density non-uniformity takes place. 

In case a fission density maximum is at the boundary of the reflec- 
tor and the core due to a reflector moderating effect U-235 displacement 
from peripheral regions of the oore to central regions results in both 
improvement of the heat generation non-uniformity and the reaotor charge 
gain. 

In a thermal heterogeneous reaotor redial flattening of the average 
energy generation in separate fuel elements may be reached by gradual mo- 
derator volume per cent rising from centre to the periphery of the oore. 
Parallel with this it is possible to uniform U-235 efflolenoy function 
through fuel elements. As the energy release non-uniformity problem is 
of importance In engineering, the variational problem of orltioal reao- 
tor mass minimization with the fission density non-uniformity ooeffloi- 
ent limiting in the oore may be formulated. 

The derivation on Euler’s equation for this problem is identical 
to one mentioned above. In deriving this equation one takes into aooount 
an Isoperlmetrio condition for the voliune non-xml form! ty ooeff lolent 
while an expression for /V may be obtained on the baela of 
the formulas In Ref /8/. As a result the oonstanoy condition of U-235 
efflolenoy function is obtained in the oore, having a special expression 
for the problem under oonsideration. 

§5. Variational problem for niazimum of breeding ratio 
(BR) in a breeder. 

Consider a breeder with given demenalons and given oore and r^jfleo- 
tor materials. We assume the reaotor fuel is U-230 and Pu-239 and there 
Pu— 239 is aooumulated. Aooordlng to the detomlnation the breeding ratio 
In a reactor throughout la; 

RD - _£• (£) c/B dftdv^ 

' “/f ” /// F(=c,e,K) Ps C'2) o'i? rfv (8) 

The relation bMween the U-238 oonoentratlon />gC7) end th« Pu-239 oon- 
oentretion />g(?) in every raeotor point r" can Ije varied v/lth the follo- 
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wln^ condition preserved; 

r J3 ( 9 ) 

irhere andys are some given fuel element parametero. What distribution 
of concentrations J> ^Cf) and /3g(r) in the reactor ie required to achie- 
ve the maximum BR without variation? 

Expressions for the variation of the baoio funotional BR and ieope- 
rlmetrlo condition <f have been taken from Ref/8/, let <7g(f) be a 

varying funotlom because/) g(f) is connected with it by the relation (9). 
Taking into aooount oonsidoratlon mentioned above tre vrrite; 
p,)Sp,(i:)dv 

Sx*rr )<^Pt(t)dv=o ^ ^ ^ 

where funotions ) and f (f;o ) are as follows: 

*[il (E ') e>* (£ V Jf « (£, £;sn')~i- v^/r£ 'je/ (e ') x, (£.e: sin') ]f » 
^P'(i.£:sJ:')<fy*(^.E,S.}d£'cla'dEdS ' (12) 

. oA'oii ^ 

ni/(E') 6 ;(E')XsiE.E:slSi')- ^ i/(E')X>(E.E: 

E*(%E,3:) dE'dK'ctEdn - 

Z (d<,tfs(£)] F'(z.E,n)FUt,E.ii)JEdst (13) 

By moans of relations (10) and ( 11 ) we found Euler's equation for the 
problem und^r consideration • 


e(=^/p,) = tlijM. 


- const 


The condition obtained of constancy U-238 efficiency function in the 
reaotor for problem under oonslderatlon can be used in numerloal calcu- 
lations to find the U-238 optimum distribution p g(r). 

Chapter III, OPTIMIZIBG OP PHYSICAL CHARACTERISTICS 
OP RADIATION SHIELDING, 

§1, The statement of problem. 

The basic problem of design and oaloulation of reaotor shielding Is 
to find suoh distribution of shielding materials (heavy and light oompo- 
nents, neutron absorber eto.), that satisfies some requirements. These 
requirements ussually consist in minimising of a shielding weight or sl- 

ses under definite oondltlons (given dose on the shielding, heat genera- 
tions eto,). 
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Thus, this problem may be formulated as a variational problem on 
optimum of some value for other values given, 

In some simple cases this problem can be resolved as a olaasical 
variational one (E.Elizard, L.Klmel ,E. Petrov, L. Pankratov, G.Lisoohkln et 
al /20/,/21/,/22/), However, in more complex cases when a large change 
of radiation spectrum occurs (for instance, while introducing strong 
absorbers for capture gamma radiation depression), processes in shield- 
ing should be described by a complex multigroup system of equations and 
using of classical methods is difficult. Classical methods apply also to 
problems, when optimum is achieved in an internal region of pemilssible 
material concentration alterations (r)^ where p (?) is concentra- 

^ q ^ ' q ^ 

tion of "q” material in point r^p ^ is its maximum value). However, in 
practical problems there is a situation, when in some regions it is not 
advantageous to put some materials, and In other regions it is profitable 
to put a material with maximum dencit.y, 

A general approach to the problem solutioi; on optimum distribution 
of shielding materials may be based on using the method of consecutive 
approximations, step "n’' of which is a transition from distribution 
p^^^(r) to P Functions should be 

chosen in such a way that : 

a) the minimized value would decrease as much as possible; 

b) V’alues F^,F2-«*Fjjj v/ould remain equal to given ones; 

o)tranBition would remain in region of permissi- 

ble values of functions y) ^(F^ . 

To choose cfp ^^^(r) it is necessary to know in what way the value 
will be changed if the density of the materials changes in different 
points ?• If the value is rather small for the definition of 

variation (5^ formulas of small perturbation theory (Refs /7/,/10/) can 
be used, where F^ is expressed by means of functions of flux and "dan- 
ger” of neutrons, which are solutions of the basic equation and the ad- 
joint equation in the preceding iteration* These formulas permit to cal- 
culate functions: 

iPAf) 

which are material "q" effloienoy in the point r relative to the value 
The knowledge of efficiency functions permits reasonable choice of va- 
luecTp'^^^C^ In every iteration. For example, in the absolute extremum 
problem can be chosen in form (fy?^^^(?)aAh^^^(r'' in suoh a way 

that If q material efficiency relative to P^ is maximum, the greater 
amount of this material is added, where is negative its amouirt is 

decreased. Of course, in order to use the small perturbation theoiy it 
is necessary to choose rather small value of (?), this oan be ensured 
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by choosing a proper value A. 

For the conctructlon of numerical algorithm of optimum solution in- 
vestigation it can be used the linear programming ideas. 

§2. The perturbation theory application to shielding 

characteristics. Shielding material efficiency functions. 


In general physical characteristics of a shielding 
depend on neutron flux, gamma quantums and distribution of shielding ma- 
terials 

The neutron flux distribution is described for the steady state by 
the equation: 

where 2CE) is a total macroscopic cross-section of Interaction; 

is the differential section defining the transition of 
neutrons from the beam to the beam (here it is considered 

an elestio and inelastic transitions and fission too); q(f , JT ,E) is 
the distribution of neutron sources. A space for which we find a soluti- 
on of equation (1) occupies the volume V, limited by surfaces S* 

The different physical characteristics of a shielding can be writ- 
ten in the form: 


where the first term defines values due to the different neutron interac- 
tion processes while passing surface, a heat generation from neutrons 
and capture gammar-quantums, a neutron dose, etc. The term c^^^(?)yde- 
flnes values due to extemal gamma-quantum sources (for instanoe, gamma- 
quanturas form active core). Here we have used the approximate method to 
consider gamma- quantum, although its description by a kinetic equation 
should not be a principal difficulty. 

The shielding weight is a functional in form: 

an adjoint with (1' 


Consider an 


1) equation: 


where is a functional derivative of with respect to * 

In case of linear dependence from for functionals of type (2): 

kB has bfeen shown In Refe /6/,/7/,/lO/,/13/, the adjoint function 
,B) describes a variation of value if one additional neutron 
has been injected In point In shielding this function defines a 

neutron "danger” relative to the neutron dose over a shielding, to gamma^ 
radiatlon dose, to a heat generation, eto. 

A knowledge of flux and "danger" functions permits to formulate a pertur- 
bation theory for different functionals. 


It oan be shown, that 

364 - 14 - 
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where and cf ^ are variations, due to a material density variation, 

sp,m. at the kinetic equation operator (l), and the right parts of 
equations (l) and (3) too. 

As it has been shown in §1 the efficiency material function re- 
lative to different functionals has foimi: 

The efficiency function defines the functional variation F, while 

. ^ K 

inserting a unit quantity of the material in point r. Values 
can be calculated using the perturbation theory formulas (a). 

Then a change in functionals at varying material q density in 

a unit volume near a point ? by a value ^(f) will be equal to: 

Thus having informat ion^about efficiency functions know what varia- 
tions in different functionals will result from some small variation of 
the material densit ies(fj3^C?)- 

§3. The minimizing procedure. 

With consecutive operations for minimizing some functional pro- 
vided other functionals P^ . . . are under definite conditions: 

it is necessary to find such valuesap^(r) in each stage that should lower 
the value of the minimized functional under the condition (6): 

Let US write p (?) in the form 

where^ is a theoretical denr^lty of material q, 

and (7) 

called further the control, will be considered an piece-cons- 
tant and equal to at interval Afj. let ho in vector form iTj, which 
has components U .. 

Let ^ P in form whe re tP ^ ^ ^ v n. r i a t i o n s of 

vector O"^, the using of which with positive weight docs not take out 
from the region I' 7) $U t and-cfV^^ are considered to be different, some 
of them can be ” barred”;, are some numbers, satisfied by limits: 

^ ( 8 ) 

Vectors correspond to variations vectors components are 

defined above by efficiencies h^. of material in a point J relative to 
the functional 

/ectoi appcilainc to a_^npor;e , measurement of v/hich Is (mi). 

In the construction of vectors h ; 'r Ihe space it is necessary to 

introduce a metric, ns different components of vectors h . have different 

n J 

measurements and different orders of values. The metric can be introduced 
for example, by means of variation of scales on the different components 
(i.e.by action of some diagonal matrix) so that average values of vector 
components h ^ should become approximately equal. 

All further constructions are based on the 1 Inear small perturb?^t ion 
theory. Now the variation problem is ir< the foil owing form, 

36^f 
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control 


To find such numberB S^j satisfying to limitations (0), that for the 




(fily 


conditions (6) and(7}should be fulfilled and with this an increment of 
the mini ml zing functional should be negative 

A number of ^oints*^^ 

at all Sqj, subordinated by the inequalities (0), have made a convex clo- 
sed polyhedron R, which should be called an attainable region, n Is a li- 
mited region of the space R 

m+i 

To definition let us consider a case, when conditions (6) have a 
form of equalitlea. Then the problem is reduced to a search of the point 
of the intersection of the ray 1^ with the direction (-1,0,0...o) with a 
boundary of region R Cin 'tbe common case to an intersection point of some 
Bpaoe, the points of which satisfy to conditions (6) with the boundary of 
the region. 

let us consider one of simple but in practice sufficiently effective 
methods to define S j. (In detail this problem has been considered in Ref. 
/ 23 /). 

The variation procedure consists of series of the elementary actions, 
in everyone of which the control varies only on one -a xy for one component q. 
A variation is chosen in such a way that as a result the minimized functi- 
onal being more lowered at a least possible disturbance of the conditions 
(6). Such variation ohosen for the conditions of minimum angle between 
the vectors h^^^ and 1^. Let pass from U to •c/’U .. As consequ- 

eno^ of this there will be a displacement of the point f to the point 
distance of which from axis 1^ is 

The second elementary variation act Is the restoration of boundaiy 
oondltions(6). In this case the control variation is ohosen from conditi- 
on of minimum angle between vectors h^^ and (0,|\ 

All vectors h j available at our disposal are considered identically. 
In that the Iteration Is completed . vectors 


Ls completed^ vectors h^^are counted again at c^nt- 


previous iteration and then a new iteration is 

begun* 

Other method of values S^^determination is possible. Itconslsts in the 
oholoe of some distance ^^from the axis (-1,0,0..,0)(lt is some boundaiy 
condition non-oon8ervanoe),whloh tends to zero in the iteration procedure. 
Then detennlned from relation: 


p{ P =p( = S^jP(pr^^- r 


The prooess is stopped, if it is not possible to select such S by means 
of which the fimotional P can be reduced with consorvanoe of condition 
( 6 ). 

The optimality criterion oan be formulated in such a way. The optimal 
oomposltlon is one (ouch values of control U^p, vectors IT of which 
form a convex cone, not containing the minimized direct ion'^vector t . 

Denoted minimization principle is applied, of course, to reactor 
problems too* 

-, 6 - 
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